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Abstract 

In two-stage robust optimization the solution to a problem is built in 
two stages: In the first stage a partial, not necessarily feasible, solution 
is exhibited. Then the adversary chooses the "worst" scenario from a 
predefined set of scenarios. In the second stage, the first-stage solution 
is extended to become feasible for the chosen scenario. The costs at 
the second stage are larger than at the first one, and the objective is to 
minimize the total cost paid in the two stages. 

We give a 2-approximation algorithm for the robust mincut problem 
and a (7+2)-approximation for the robust shortest path problem, where 
7 is the approximation ratio for the Steiner tree. This improves the 
factors 1 + V2 and 2(7 4- 2) from [Golovin, Goyal and Ravi. Pay today 
for a rainy day: Improved approximation algorithms for demand-robust 
min-cut and shortest path problems. STACS 2006]. In addition, our 
solution for robust shortest path is simpler and more efficient than the 
earlier ones; this is achieved by a more direct algorithm and analysis, 
not using some of the standard demand-robust optimization techniques. 
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1 Introduction 

The general setting in a two-stage optimization problem is as follows: There 
is a set of demands (aka scenarios), one of which has to be satisfied tomorrow. 
It is not until tomorrow that it is revealed which demand must be satisfied. 
A demand is satisfied by buying a set of resources. Thus, one possibility to 
satisfy the tomorrow's demand is to wait until tomorrow, know the scenario, 
and buy a corresponding set of resources. However, the resources are cheaper 
today than tomorrow, by an "inflation" factor A > 1. Hence it makes sense 
to buy some resources already today, i.e., at the first stage, even without 
knowing the tomorrow's scenario. (Say, if A = 00, the resources bought 
today should better satisfy all demands.) Then tomorrow, upon revealing 
the requested demand, only some additional, second-stage, resources have to 
be bought at the higher price. 



The tomorrow's demand is chosen by an adversary. The adversary knows 
the resources bought at the first stage. He also knows the algorithm that 
you will use for buying second-stage resources. The adversary chooses the 
scenario so that your second-stage cost is as large possible (the adversary 
is omnipotent, and can solve an NP-hard problem for that, if necessary). 
Your objective is thus to minimize the maximum, worst-case cost paid in the 
two stages. Because of such hedging against the worst demand, this type of 
two-stage robust optimization is called demand-robust. 

Related work 

In stochastic optimization (see, e.g., [8, 9j) the objective is to minimize the 
expected cost paid over the two stages. Universal approximations [U [TO] 
may in a sense be viewed as one-stage robust solutions. The demand-robust 
optimization as studied in this paper was introduced by Dhamdhere et al. 
in PJ. Several techniques have proved to be viable in the field: 
Minimal feasible solutions. Dhamdhere et al. j2| showed that there always 
exists an approximate first-stage solution which is a minimal feasible solution 
for a subset of scenarios. Restricting oneself to such solutions makes one 
loose at most a factor of 2 in comparison with the unrestricted case. Since 
the pioneering paper [2], the minimal-solution idea has been extensively used 
in the design of approximation algorithms for two-stage robust optimization 
problems. 

LP rounding. IP formulations of optimization problems often extend 
directly to stochastic and demand- robust versions; rounding the LP relaxation 
solution is one way to obtain an approximation. 

Thresholded a- approximations. A common approach to demand-robust 
optimization is as follows: Suppose you are shooting for an a-approximation. 
Guess the second-stage cost C* n of the optimal solution (often the number of 
relevant CjjS is small; if worse comes to worst, go through "all possible" C^s 
approximately with repeated doubling - or more precisely, with repeated 
(1 + e)-ing). In the first stage, satisfy all high-cost demands - those each 
of which is more expensive than aC* u to satisfy. Then in the second stage 
you are guaranteed to pay at most ctC* II - which is within factor a of 
optimal second-stage cost. Finally, argue that your first-stage solution is 
also within a times the first-stage cost of the optimum - for the overall 
approximation guarantee of a. A very general treatment of the thresholded 
covering algorithms is presented in a recent paper [7]. 

Our contributions 

In Section [2] we present a thresholded 2-approximation for the robust mincut 
problem. This improves the (also thresholded) (1 + \/2)-approximation from 
[3j . The improved approximation guarantee is due to a refined analysis using, 
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similarly to [4J, laminarity of mincuts (the Gomory-Hu mincuts tree). 

In Section [3] we give a (7 + 2)-approximation algorithm for the robust 
shortest path problem, where 7 is the Steiner tree approximation ratio. 
This improves the 2(7 + 2)-approximation from [4J (the techniques in [6] 
potentially imply a 4.25-approximation). The algorithm and its analysis are 
very simple. In particular, unlike [3] we do not restrict ourselves to minimal 
feasible solutions and do not use the thresholding. Avoiding the guessing of 
the second-stage cost of the optimum makes our algorithm more efficient (by 
at least a linear factor) than that of [4J. 

2 Demand-robust mincut 

In the demand-robust mincut problem the input is a (positively) weighted 
undirected graph G = (V, E, C) with C representing the capacities of edges, 
the root vertex r E V, and a set T C V \ r of terminals. For a terminal t £ T 
let m(i) denote the minimum r-t cut (if the mincut is not unique, take m(i) 
to be the cut that cuts out from r a maximal set of vertices); similarly, for a 
set S C T of terminals, m(5) is the minimum r-S cut. We use C(t),C(S) 
to denote the capacities C(m(i)), C(m(S)) of the mincuts m(t), m(S'). For a 
subset E 1 C E 1 of edges let m (t) be the minimum r-t cut in G with weights 
of edges in E' set to 0; let C^i(t) denote the capacity of the mincut m#/(t). 

A feasible solution to the robust mincut problem is an arbitrary set 
Ei C E of edges. The cost of the solution is 

C(E[) + A • maxCWt) 

where A is the inflation factor. 

The edges Ej of the solution are called the first-stage edges and the 
cost C(Ej) — first-stage cost; the edges m^ J (t) are the second-stage edges 
for terminal t and the cost max^gT CEj(t) is the second-stage cost. The 
objective is to find Ej minimizing the two-stage cost (with the second-stage 
cost inflated by A). 

Notation For a set P C V of vertices let dP denote the boundary of P — 
the set of edges that have exactly one endpoint in P. We use Ej to denote 
the optimal solution. 

2.1 Mincuts laminarity 

Let G* be G with the capacities of edges in EJ set to 0. For a terminal 
t £ T, let Q* t C V \ r denote the t-side of the cut m e* (t) — the vertices 
reachable from t after the edges Ej and iriE*(t) are removed (the asterisk 
emphasizes that Q* is the t-side of the mincut in G* , not in the original G) . 
It is known (e.g., can be seen from the Gomory-Hu tree Section 8.6]) 
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that these i-sides do not properly intersect — Vu, v £ T either Q* u n Q* v = 
or Q* u C Q* or Q* C Q*. In other words, for any subset 5 C T of terminals 
the t-sides of the terminals in S form a laminar family J 7 ^ = {Qjf : t £ 5} of 
sets. 

Let Fg C J 7 ^ be the basic (inclusion-maximal) sets in the family J 7 ^.; 
assume that all sets in Fg are unique (note that in principle we could have 
Qu = Qv = Q* f° r u,v £ S,u ^ v, with Q* not being a proper subset of any 
other set in F s — in this case only one of Q* , Q* v is included in Fg). Call 
the terminals B s = {b £ S : Ql £ Fg} the basic terminals of S. 

2.2 Thresholded a-approximation 

The thresholded covering paradigm applied to the robust mincut problem 
works as follows: Imagine that we know the cost Cjj = maxjg^ Ce* (t) that 
the optimum pays at the second stage. To obtain an a-approximate solution, 
cut out the set U = {t £ T : C(t) > aCjj} of "expensive" terminals in the 
first stage. That is, the output of the algorithm is m(U). 

Assuming the terminals in T = (t\, . . . , t\j>\) are ordered in non-increasing 
order of mincut (C(tj) > C(t; + i)), for any C* n we have U = {t±,t2, ■ . ■ ,tj} 
for some j = j{C* n ) £ {0, 1, 2, ... , \T\}. Hence there are only \T\ + 1 different 
possible sets U for all possible C* n > — so try all the possibilities and 
choose the best; this way the non-determinism in C is reduced just to the 
non-determinism w.r.t. j. 

For the approximation ratio analysis assume that the algorithm is run 
with the "right" guess of U corresponding to the right choice of Cf T . By 
definition of U, the thresholded algorithm pays at most aC* H in the second 
stage. The tricky part is to bound the cost, C(U), of the first stage. 

The analysis of |4j Golovin, Goyal and Ravi [1] use the following estimate 
of the first-stage cost of the thresholded a-approximation: 

C(U) < C{ [J dQt) < (1) 

beB* v 

< C{E}) + ^— C(dQtnE*j) 
a — 1 *■ — ' 

beBi T 

To bound the last sum, [1] cleverly use the fact that due to pairwise- 
disjointness of the basic sets i 7 ^, every edge e G Et appears at most twice in 
the sum; thus C(U) < (1 + ^j)C(Ej), and the overall approximation ratio 
of the thresholded a-approximation algorithm is max(l + a), minimized 
by a = 1 + v2 — the final approximation ratio of [1] . 
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Using a = 2 is enough 

Our algorithm is just the thresholded 2-approximation, i.e., the output of 
our algorithm is the minimum r-U cut where U = {t £ T : C{t) > 2C|j}; as 
usual, for the analysis we assume that Cjj (or, equivalently, U) was guessed 
correctly. The second-stage cost of our solution is at most 2C* n . In what 
follows we prove the bound C(U) < 2C{E* I ) on the fist-stage cost of our 
algorithm. 

To show C(U) < 2C(Ej), instead of a correct but too generous bound 
([!]) of |4J on C(U), we use a tighter estimate 

C(U) < C Id |J Ql 
\ beB* v 

The correctness of the estimate follows from the same argument as in 
every terminal of U belongs to at least one of the sets Ql , and none of the 
sets Ql contains r; thus the boundary of the union is an r-U cut. 
To prove 



C id |J Ql) < 2C(E*j 



we argue that 



C id |J Qt) < 2C A* |J QJ (2) 

where A*(P) = {{P x P) n E*j) L> (dP D E*j) denotes the edges from E} 
that have at least one endpoint in a set P C V of vertices; clearly, the 
right-hand-side of Q is at most 2C(E|). 

Number the terminals in arbitrarily: By = (61,62,...). For k = 
0, 1, . . . , \Bfj\ define B* k = Q* bi U Q^ U • • ■ U Q* bk . The inequality g follows 
from the next lemma: 

Lemma 2.1. VA; = 0, 1, . . . , \Bfr\, C (dB* k ) < 2C (A* (5*)) 

Proo/. By induction on k. The base is trivial: = (7(50) < 2C(A*(0)) = 0. 

Let X* = d(B*_ v Q* bk )n E*j,X = d(B* k _ v Q* bk )\ E*,Y* = (dQ* bk \ (X* U 
X)) n E*j, Y = (dQ* bk \ (X* U X)) \ E} (Fig. [1]. 

Because bk belongs to the set of high-cost terminals U, the optimal 
solution Ej must "help" 6^ by at least half (using the terminology from [7], 
b k is "low"): 

C(X*) + C{Y*) >C{X) + C{Y) (3) 

Indeed, since dQ* h is an r-b k cut, C(b k ) < C{dQ* bk ) = C{X*) + C(Y*) + 
C(X) + C(Y), and since 6 fc G U, C{b k ) > 2C* n >2C E *(b k ) = 2(C(X) + 
C(Y)), from where ^ follows. 
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Figure 1: E* is bold. B* k = B*_ x U Q* 



Using C(X) > we rewrite ^ as 

C(Y*) > ^1 + C(Y)-C(X*) + C(X) > 

> C(Y*) + C(Y)-C(X*)-C(X) (4) 

We have (see Fig [I]): 

C(A*(B k )) > C(A*(Bl_ 1 )) + COT) (5) 
C{dB* k ) = CidBl^) + C(Y*) + C(Y) - C(X*) - C{X) (6) 
By the inductive hypothesis, 

C(A*(jBfc_i)) > C{dB * k ~ l) ( 7 ) 
Putting @, ([5j), together we obtain 

c(A*(B k )) > cidB ^ | g(r*) + c(y)-c(x*)-c(x) = c(dB* k ) 

□ 

Overall, we have that the first-stage cost of our solution is at most 2C(Ej), 
and the second-stage cost is at most 2C* n : 

Theorem 2.2. There is a polynomial-time algorithm which gives a 2- 
approximation for the robust mincut problem. 



3 Demand-robust shortest path 

The input to the demand-robust shortest path is the same as to the demand- 
robust mincut problem: graph G = (V, E, w) with w representing the lengths 
of edges, root vertex r £ V, and a set T C V \ r of terminals. A solution is a 
set Ei C E of edges. The cost of the solution is 

w(Ei) + A • maxw(SPE I (t)) 

where SP^; / (t) is the shortest r-t path in G with weights of edges in Ei set 
to 0. The objective is to find Ej minimizing the two-stage cost. 
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3.1 Algorithm 



Our solution is a Steiner tree on a subset S C T of A terminals. The set S is 
built incrementally, starting from r, and repeatedly adding a farthest (with 
ties broken arbitrarily) terminal, until gathering A of them: 

1 S <- {r} 

2 while IS"! < min(A, \T\) 

3 do S 1 S U arg max SP(ii, S) > Add farthest terminal 

4 return £7 approximate Steiner tree on S 1 



3.2 Approximation ratio 

We now analyze the approximation guarantee of the algorithm. Let / be the 
distance from S to the terminal added in the last iteration of the while loop 
(line [3]). Because \S\ was growing from iteration to iteration, we have that 
at any iteration the distance from S to the farthest terminal was at least 
/. Thus, the distance between any two vertices in S is at least /. Hence, 
St(S') > \S\f/2 where 51(5) is the weight of the minimum Steiner tree on 
S (to see this, go twice around the tree — you traveled 2St(5), spending at 
least / traveling between any two vertices in S). That is, 

f < 2St(5) 
J " \S\ 

Because we always add farthest terminal to S, at the completion of the 
algorithm the distance from any terminal to S is at most /. Hence, 

w(E H ) < f 

where En are the edges that we buy at the second stage. 

Let Ej be the optimal solution. Let Ejj(t) be the edges that the optimal 
solution buys at the second stage if the demand is t G T. Let E* n be the 
edges that the optimum buys in the worst case: w{E* n ) = maxj w(Ejj(t)). 
Then Ej U te g\ r Ejj(t) is a connected graph that spans S. Hence 

w(E* r ) + ™(Eh(t)) > St(S) 

teS\r 

and 

w(E*j) + (\S\ - l)w(Eh) > St(S) 

We consider the cases A > \T\ and A < \T\ separately: 
If A > \T\, then S = T U r, and 

opt = w (Ef) + Xw (E* n ) > w{Ej) + \T\w(E* n ) = 
= W (EI) + (\S\ - l)w(E* n ) > St(S) > w{Et)h = apx/ 7 
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where opt is the optimal cost, apx is what we pay, and 7 is the approximation 
factor for the Steiner tree. That is, apx < 7 opt. 



because A < \S\ < A + 1. 

Theorem 3.1. If for some class of graphs there is a ^-approximation for 
Steiner tree, then for that class of graphs there is a (7 + 2) -approximation 
for robust shortest path. 

For general graphs, the best 7 = 1.39 is due to Byrka et al. jTj: 

Corollary 3.2. There is a polynomial-time algorithm which gives a 3.39- 
approximation for the robust shortest path problem. 

3.3 Running time 

As far as the efficiency of approximating the robust shortest path is concerned, 
the best bound that can be given on the running time of the algorithm of |4J 
is 0(|T||V||2£|). The (multiplicative) overhead of 0(T||V|) is due to guessing 
|T||1/| possible values for the second-stage cost of the optimal solution. Then 
for each guess the algorithm of [3j builds an approximate Steiner tree, which 
must take time. 

Because our algorithm avoids the guessing, we can achieve the running 
time of 0(min(A, iTDd^l + |V| log |V|)) at the expense of increasing the 
approximation ratio to 4. For that, in line [4] we use the 0(|-E| + \V\ log | V|)- 
time 2-approximation algorithm of Mehlhorn [12] or Floren |3j . Then our 
algorithm's running time is dominated by finding the farthest terminals in 
line [3 

Corollary 3.3. There is an 0(min(A, |T|)(|.E| + \V\ log \V\))-time algorithm 
which gives a ^-approximation for the robust shortest path problem. 
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